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Recent developments in many-body quantum chaos have raised the is- 
Q>^ ' sue of correlations between different families of levels in the spectra of 

, random fermionic systems. It seems that rotational invariance is sufficient 

■ to force an otherwise random system to have statistically ordered spectra 

I where the ground state can have zero spin or the maximum spin allowed a 

i huge portion of the time. We have tackled this question and have charac- 

terized the role of "geometric chaoticity" and pairing in random systems. 
We have shown how the geometry gives rise to such a distribution of ground 
• state spins and examined the coexistence of these regularities with chaotic 

wave functions. 



^3 . PACS numbers: 24.60.Lz, 21.60.Cs, 05.30.-d 



1. Introduction 



5h I The role of chaotic dynamics in nuclear physics has been extensively 

^ studied. Geometric chaoticity, or the assumption of pseudorandom coupling 

of angular momentum, was used to derive the early level density formulae. 
Random matrix theory clarified the role of global symmetries and lead to 
new methods of statistical analysis of nuclear spectra, while work with shell 
model wave functions gave a deeper understanding on the thermodynam- 
ics of the nucleus. The question of correlations between different classes 
of states described by the same hamiltonian has been neglected. Recent 
investigations [1, 2] showed that the spectra of random hamiltonians which 
respected global symmetries exhibited strong correlations between the dif- 
ferent classes. This phenomenon was subsequently observed in a variety of 
systems but recently a reasonable physical explaination was proposed [3]. 
This work addresses this important issue and explains the correlations in 
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terms of geometric chaoticity. The preference of the systems under study 
to have spin zero or maximum spin ground states and parabohc yrast hnes 
is shown to have its origins in the geometry of the system, as opposed to 
the specific details of the interaction, i.e. it is a direct result of the com- 
plex coupling of many single-particle angular momenta required to build a 
many-body state with the appropriate symmetry. 

In order to address the regularities in the spectra of random systems 
we choose to study the simplest possible system, that of N fermions on a 
single energy level of spin j which has a 2j -|- 1 degeneracy. The regularities 
are robust with respect to the choice of both system and ensemble. Within 
this space, the general two-fermion rotationally invariant Hamiltonian can 
be written as 

H = J2VlPLPla, (1) 

LA 

where the operators for a pair of fermions with total angular momentum L 
and projection A are defined as 

~ ^mn^mQn) -PlA = — 1= C'mn^n^m; (2) 

V ^ mn V ^ mn 

Um annihilates a particle with single particle angular momentum j and pro- 
jection m, are the Clebsch-Gordan coefficients {LA\jm;jn) , and the 
two-particle states \2;LA) = -P^aI^) properly normalized. Because of 
Fermi statistics, only even values of L are allowed in the single- j space. 
There arc j + ^ interaction parameters Vl and these numbers define our 
ensemble. The ensemble we chose was a uniform distribution of Vl between 
-1 and 1, which sets the energy scale. We analyzed many realizations of 
our ensemble for various values of particle number, A'^, and single-particle 
angular momentum j. Spectra were obtained by direct diagonalization of 
the hamiltonian matrix for a given set of Vl but as the dimensionality got 
too large the shell model code OXBASH was used to get the spectra and 
wave functions. 

In the A'^ = 4 and 6 particle cases there was a huge fraction of systems 
with ground state spin Jq = and Jq = Jmax where Jmax is the maximum 
allowed value of angular momentum. Wc will call these fractions /o and 
fjmax respectively, see Fig. 1. The predominance of these values of angular 
momentum occurs in spite of the fact that they account for a small fraction 
of the total states. In the N = 5 particle systems there was a large fraction 
fj of ground states with total spin J = j which is consistent with the picture 
of a valence spin j particle on a spin core. 

The mean yrast line for all spectra with a spin ground state was a 
parabolic function of J with a strong odd-even effect, see Fig. 2. For the 
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Fig. 1. The fraction of ground states with J = for A'' = 4 and different j. 

Ensemble results; solid line. Statistical theory; dotted line, panel (a). Upper limit 
for from the statistical theory and analysis of the Jmax region; dotted line, 

panel (b). 



maximum spin ground state spectra the parabola is inverted. These smooth 
yrast lines show a strong correlation between levels with different quantum 
numbers. These correlations coexist with the familiar ones of random matrix 
theory. The distribution P(s) of nearest neighbor level spacings within a 
set of levels with the same angular momentum was in full agreement with 
the Wigner distribution. 



2. The Statistical Model 

The system was treated with an approach based on equilibrium statis- 
tical mechanics which leads to the Fermi-Dirac distribution 



exp(7m — /x) + 1 

under the constraints 

N = Y,nm, M = J2mnm. (4) 

m m 

The quantities //(iV, M), the chemical potential, and j{N, M), the cranking 
frequency, are the Lagrange multipliers associated with the two constraints. 
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Fig. 2. The mean lowest energies vs J for the N = 6, J = 11/2 ensemble for spectra 
with (a) Jo = and (b) Jq = Jmax- Separate values of the moment of inertia I 
were extracted from the odd and even J points.. 



The equilibrium energy for an N particle-state with total projection M is 

2 



r<LA 

^12 



{nin2). 



(5) 



{H)= E 

L,A,1,2 

Treating 7 as an expansion parameter, expanding rim, performing the sum- 
mations over the various products of single particle projections and Clebsch- 
Gordan coefficients, and identifying AI with J wc come to 

■4 



{H)n,j = ^{2L + 1)Vl ho{L) + h2{L)J^ + hi{L)j' 



(6) 



where, with ft = 2j + 1, 

/.„(L) = ^, ML) = ^(L=-2f). 
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Fig. 3. The average values of the random parameters Vl for the ground state spin 
Jo = Jmaxi N = 4, and j = 17/2 as a function of L shown with the statistical 
errors (solid line), and the results of the statistical model for /i2 < (dashed line). 



The regularities of the ensemble are insensitive to the choice of distribu- 
tion of interaction parameters Vl as long as the choice is symmetric about 
zero. The distribution of the sum in Eq. (6) can be calculated assuming a 
Gaussian distribution for each Vl- Our model captures the upper limit for 
/o in the systems of various j, Fig. 1. The staggering of /o in the numerical 
results for N = 4, with its period of 3 in j, can be understood by observing 
the behaviour of the multiplicity of the spin states as j changes. It is 
incremented by 1 for every time j is incremented by 3. The theory overes- 
timates the value of fjmax but the expansion of the 7 parameter was valid 
only for small values of M and thus J. 

We can use (6) to extract the mean values (Vl) of those sets of interaction 
parameters that generate spectra with certain Jq. In Fig. 3 we compared the 
theoretical points for {Vl) vs. L for AT = 6 systems with the actual mean 
values of random parameters Vl which led to Jq = Jmax- The resulting 
curve is in total agreement with our statistical theory. 

The equilibrium energies from Eq. (6) are highly correlated with the nu- 
merical values, differing by a constant negative shift. When plotted against 
each other the points fell along a straight line with a slope of 0.81 it 0.02 
and 0.81 it 0.03 for those systems with Jo = and 30, respectively, for all 
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Fig. 4. Numerical values of the lowest J = 2 (a), and J = 27 (b), energies vs. 
predictions of the statistical model, eq. (6), for the Jo = spectra of the N = 6, 
j = 15/2 ensemble. 
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Fig. 5. The constant shift in the theoretical equilribium energies of spin J vs. J 
for the N = 6 j = 15/2 ensemble for the two values of ground state spin. 



values of J, see Fig. 4. There was an odd-even staggering in the value of 
the negative shift and its magnitude was smaller for Jq = 30, Fig. 5. 



mulhall printed on February 8, 2008 



7 



0.5 




CD 0.6 
> 




-8 -6 -4 -2 2 4 6 8 



m (single particle projection) 

Fig. 6. The average occupation numbers for the 1417 out of 2000 (70%) spectra 
with Jq = Q for the N — 6, j — 15/2 ensemble. 



The assumption of statistical equilribium implies that the off-diagonal 
matrix elements of am and aj„ and their products make an incoherent con- 
tribution to (6). We also assumed that the occupation numbers are uncor- 
related. Actually using the numerical results for the correlated occupation 
numbers in (5) resulted in a slope close to unity in the plots of Fig. (4). In 
Fig. 6 we see the average values of rim for the lowest J = 2, 4, 6, 8 energy 
levels of those N = 6, j = 15/2 spectra with Jq = 0. They oscillate about 
the statistical values, exhibiting a shell structure that the theory smeared 
out through the approximation of uncorrelated occupancies. 

In spite of the regularities of the ground state spin, their origin in random 
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Fig. 7. The distribution of overlaps, eq. (8), of ground states with Jo = with the 
fully paired (seniority zero) state, solid histograms, for the ensemble with all Vl 
random, panel (b), and for the ensemble with random Vl^o and regular pairing, 
Vl=o = — 1, panel (a). Dashed histograms show the predictions for chaotic wave 
functions of dimension d = 3, panel (b), and d = 2, panel (a). 



spin coupling leads us to expect the corresponding wave functions to be 
chaotic. As a reference point, we look at the distribution of overlaps, x, 
between the ground state wave functions of our ensemble and the spin zero 
seniority zero wavefuntion which appears as the ground state for the regular 
pairing interaction, Vq = — 1, Vl^o = Oj 

x = |(J = 0,g.s.|0, p)|2. (8) 

The histogram of the distribution V{x) is presented in Fig. 7(b) for the case 
of iV = 6 particles on the j = 11/2 level where there are 3 states of spin 
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J = 0. In the chaotic hmit this would lead to 

Vd=z{x) = TT^, < X < 1. (9) 

Clearly, the observed overlap agrees qualitatively with the chaotic limit. 

Another situation is shown in Fig. 7(a) where the data, a solid his- 
togram, are taken for the ensemble which contains regular attractive pairing, 
Vq = — 1, plus uniformly random interactions in all channels with L > 0. 
The presence of pairing generates a significant probability for the paired 
ground state, a peak at x = 1. However there exists an "antipaired" state 
with J = but a vanishingly small probability to be a ground state so 
that the space of the candidates for the ground state position is effectively 
two-dimensional, and we have 

rd=2ix) = , }^ , , < X < 1, (10) 

1T\/X[1 — X) 

which again agrees qualitatively with the data, Fig. 7(a). 



3. Conclusion 

The regularities in the spectra of angular momentum conserving random 
two-body interactions in systems of N fermions on a single-j level are the 
result of geometric chaoticity, the inherent complexity of the A^-body spin- J 
states from the coupling of spin-j particles. Treating the system as a Fermi 
gas in statistical equilibrium explained the trends both in the ground state 
spins and the yrast line. The role of pairing is minimal. The shell structure 
evident in the single particle occupancies sheds light on the limitations of the 
theory. This work was supported by the NSF grants 96-05207 and 00-70911. 
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